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ABSTRACT 

We investigate the effect of a new triple-tr reaction rate from Ogata et all (|2009) on helium ignition conditions 
on accreting neutron stars and on the properties of the subsequent type I X-ray burst. We find that the new rate 
leads to significantly lower ignition column density for accreting neutron stars at low accretion rates. We 
compare the results of our ignition models for a pure helium accretor to observations of bursts in ultra-compact 
X-ray binary (UCXBs), which are believed t o have nearly pure helium donors. For m > 0.001 mEdcU the new 
triple-LT reaction rate from lOgata et al.l (12009) predicts a maximum helium ignition column of ~ 3 x 10 9 g cm -2 , 
corresponding to a burst energy of ~ 4 x 10 40 ergs. For m ~ 0.01 mndd at which intermediate long bursts occur, 
the predicted burst energies are at least a factor of 10 too low to explain th e observed energies of such bursts in 
UCXBs. This finding add s to the doubts cast on th e triple-a reaction rate of Ogata et al. ( 2009) by the low-mass 
stellar evolution results of iDotter & Paxtonl (120091) . 
Subject headings: stars: neutron — X-rays: binaries — X-rays: bursts 



1. INTRODUCTION 

The triple-alpha (3a) reaction, 3 4 He — > 12 C + y, is of 
fundamental importance for stellar evolution and explosive 
phenomena involving neutron stars and white dwarfs. The 
3a reaction rate is generally well understood at temperatures 
higher than ~ 10 K where the reaction is resonant. At 
low temperatures, the rate is dominated by non-resonant re- 
actions and uncertainties exist. Recently, lOgata et aD (120091) 
(OKK09) studied the 3a reaction rate at low temperatures 
(10 7 K < T < 10 9 K) by directly solving the three-body 
Schrodinger equation. They found tremendous enhancement 
of the 3a rate compared to previo us results, including the most 
recent rate o flFynbo et all 00051) (Fy05fl Compared to Fy05, 
the OKK09 result is a factor of 10 26 greater at 10 7 K, 10 6 
times greater at 10 8 K, and becomes equal at T > 2.5 x 10 8 K. 

The increased 3a rate allows helium ignition to occur in 
a lower temperature and density environment. This has the 
potential to strongly affect a broad range of astrophysical 
phenomena, including, but not necessarily limited to, helium 
burni ng in stars and on accre ting white dwarfs and neutron 
stars. IDotter & Paxtonl ((2009) have shown that the new rate 
affects significantly the evolutionary track of low-mass stars 
of 1 and 1 .5 M Q , resulting in a shortened or even almost dis- 
appearing red giant phase; a result that i s inco mpatible with 
observations. ISaruwatari & Hashimoto! (120101) have studied 
the effect of the new rate on the ignition of helium accreting 
C-O white dwarfs and found that off-center helium detonation 
rather than central carbon ignition may dominate the mecha- 
nism of type la supernovae, contrary to previous results. 

In this work, we investigate the effect of the new 3a reac- 
tion rate on the helium ignition conditions on accreting neu- 
tron stars and on the properties of the subsequent type I X-ray 
burst (XRB). We find that the He ignition layer is rather sen- 
sitive to the choice of the 3a reaction rate. The new rate leads 
to significantly lower ignition column density at low mass ac- 



cretion rates. We present the results of our ignition models 
for pure helium accretors and compare them directly to obser- 
vations of intermediate long bursts from ultra-compact X-ray 
binaries (UCXBs), which are believed to have nearly pure he- 
lium donors. We find that the new 3a reaction rate predicts 
burst energies that are too low to accommodate the observed 
energy of ~ 10 41 ergs of intermediate long bursts. 

2. HELIUM IGNITION CONDITIONS 
The energy generation rate of the 3a reaction is given by 



S3 a =N A r3aQ 
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where Y is the mass fraction of He and Q->, a 7.275 MeV is 
the energy release from the reaction. ri a is the 3a reaction rate 
per unit volume and (aaa) is the temperature-dependent part 
of the rate (in units of cm 6 s _1 g~ 2 ). / sc is the electron screen- 
ing enhancement factor which is most important in the high 
density, low temperature regime. We calculate / sc from the 
screening potentials by using / sc = e\p(h aa + h a B e ). Here, h aa 
and /j ff Be correspond to the screening potentials for each step 
of the 3a reaction, 4 He+ 4 He <-> 8 Be and 4 He+ 8 Be — > I2 C+y, 
respectively. For two species i and j that create a com- 
pound nucleus ({Z c , A c j = {(Z,- + Zj), (A, + Aj)}) in a reaction, 
h u = /(Tj) + f(Tj) - /(T c )0 Here, F, = (Z ie ) 2 KafaT) is the 
Coulomb coupling parameter for species i with a, being the 
ion sphere radius. /(F,) is the Coulomb free energy per ion 
for a one-compone nt plasma of species i for wh ich we adopt 
the analytical fit of iPotekhin & Chabrierl ( 120001) that matches 
with Monte Carlo simulations. To examine how the He igni- 
tion curve depends on the choice of screening enhancement 
factor, we also include in our calculations a different form al- 
ism of screening enhancement taken from lltoh et al.l (119901) . 
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1 The Fy05 rate, currently adopted by the JINA REACLIB database 
I http://groups.nscl.msu.edu/jina/reaclib/db/ I, is within a f actor of 10 of the 
previous standard rates from Cauahlan & Fowler 11988) and the standard 
rates in the NACRE compilation <Angulo et all 19991) . 



2 Here, we apply the linear mixing rule (LMR) for a multi-component 
plasma. Recent studies have shown th at the LMR is n ot accurate at the Deby e 
limit at T -» 1 IPotekhin et"aT]|2009bD and at T » 1 IPotekhin et"aT]|2009d) . 
The correction to the LMR is small in our interested range of T and is ignored 
here. 
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Fig. 1 . — Ignition curve for a pure He neutron s tar accretor. Two 3a reac- 
tion rates are adop ted: Fy05 IFvnbo et ai]|2005l) (solid lines) and OKK09 
IQgata et"afl 120091) (dashed lines). We multiply bo th rates by an elec- 
tron sc reening enhancement factor, provided e ither by|Potekhiii & ChabrieJ 
<200Q|) (thick lines) or by lltoh et alj (19901) (thin lines). The ignition 
curve with the analytic al expression for the screened 3a reaction rate from 
Fushiki & Lamb 1 1987) (dotted line) is also shown for comparison. 

When matter is accumulated on the surface of a neutron 
star, it is heated by gravitational energy release, nuclear reac- 
tions, and energy flux emanating from the deep neutron star 
crust. The temperature eventually becomes sufficiently high 
to lead to explosive burning. We determine the He ignition 
condition that is fulfilled when heating by nuclear burning be- 
comes faster than radiative cooling, i.e., 



ds 3a 
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cool 
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where e coo i = pKT/y 2 is a one-zone approximation of the 
local cooling rate. The temperature derivatives are taken at 
constant pressure P = gy. Here g is the gravitational acceler- 
ation on the neutron star surface and y is the column density, 
i.e., the accumulated mass per unit area. K is the thermal 
conductivity. We use the sa me equation of st ate and thermal 
conductivity as described in lPeng et al.l (12007b . 

Throughout this letter, we assume a canonical neutron star 
mass of M = 1.4 M Q and a radius of R = 10 km. The surface 
gravitational redshift is z = [1 - 2GM/(Rc 2 )]- 1/2 - 1 = 0.31 
and the gravitational acceleration is g — (GM/R 2 )(l + z) — 
2.43 X 10 14 cms- 2 . 

To obtain the ignition column density for a given ignition 
temperature, we increase the column density until the igni- 
tion condition (eq. (3)) is reached. In Figure Q] we compare 
the He ignition curves obtained with the 3a reaction rates 
of Fy05 (solid line) and OKK09 (dashed line). We multi- 
ply both rates by a n electron screening e nhanc ement factor, 
provi ded either by iPotekhin & Chabrierl (120001) (thick lines) 
or by lltoh et all ( 119901) (thin lines). The results of the two 
screening formalisms agree quite well. Overall, the OKK09 
rate leads to significantly lower ignition column density for 
the same ignition temperature. This effect is most promi- 
nent at temperatures below ~ 2 x 10 8 K and, hence, for neu- 
tron stars accreting at low rates. Also shown in Fig. Q] is 
the He ignition curve resulting from an analytical expression 



for the screened 3a reaction rate given by iFushiki & Lambl 
d 19871) (FL87, dotted line). This rate has found broad ap- 
plic ation in the modeling of X-ray burst ignition conditions 
(e.g.|Cummindl2003MCumming et aT]|2006t iNaravan & Hevll 
120031: ICooper & Naravanl2005l) . The~FL87 results agree quite 
well with those obtained with the Fy05 prescription. Signif- 
icant differences appear only at low temperatures and high 
densities where both the uncorrected 3a rates and electron 
screening corrections are different. 

3. IGNITION MODELS FOR ACCRETING NEUTRON STARS 

Burst ignition and the properties of the subsequent burst de- 
pend on the accretion rate, the composition of the accreted 
matter, and on the interior properties of the accreting neutron 
star. The steady-state thermal structure of the accreted enve- 
lope is gover ned by the temperature and flux equations (see 
IBrownll2000l) . 



dT _ F 
dy pK ' 
dF (dT T 

— = Cpm — "ad | - £nucl , 

dy \dy y 



(4) 
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where Cp is the specific heat at constant pressure and V a( j = 
(dlnT /dlnP)s . m is the local mass accretion rate, com- 
monly normalized by the local Eddington mass accretion rate 
wzEdd = 2m p c/[(l + X)o- TH R] with X being the mass fraction 
of hydrogen and <x XH being the Thomson scattering cross sec- 
tion. For a pure He accretor, ;«Edd = 1.5 X 10 5 g ctrT 2 s -1 . 
£nuci is the energy generation rate from nuclear burning. Here, 
we are considering only pure He accretion. The energy gen- 
eration rate during the accumulation phase alone is therefore 
governed by the 3a reaction, e nuc i = £3^. We adopt and com- 
pare the 3a rates of Fy05, OKK09 and FL87. 

We integrate from the surface (y = 10 3 g ctrT 2 ) to the outer 
crust of the neutron stai0. At the surface, we use a radiative- 
zero condition, F — 4/(3/fy) [cr R r 4 - F acc J, where k « 0.2 (1 + 

X)cni 2 g is the opacity contributed by electron scattering, 
cr R is the Stefan-Boltzmann constant, and F acc = GMih/(2R) 
is the accretion flux. We employ a fuel/ash burning scheme 
and switch the composition from fuel to ash when the ignition 
condition (eq. Q) is reached. The composition of the ash is 
set to be 56 Fe, assuming complete burning. 

At the inner boundary, we set the flux F = Fb = Qbth, 
where Qi, is the energy transported outward from the deep 
crustal heating via elec tron capture, neutron emission, and 
pycno nuclear reactions (ISatolll979t lHaensel & Zdunikll 19901 
120031) . The total energy release from the crust is set to be 
Gcmst = 1.4 MeV/nucleon, but not all of this will be trans- 
ported to the surface. Some fraction of the heat will diffuse 
to the core, provided efficient core cooling by neutrinos. The 
exact value of the outward energy release Qb is sensitive to 
the core thermal structure which we do not include in our 
calculations. Several self-consistent ignition models that in- 
tegrated deeper into the core and matched th e core flux to 
the neutrino luminosity have been carried out dBrownl 120041 : 
ICumming etal]|2006t ICooper & Naravanll2005l) . These mod- 
els predict Qb 0.1 MeV/nucleon at m > 0.1 WEdd an d in- 
creasing at lower accretion rates. However, we note that the 

3 The ignition point is not sensitive to the choice of the column density at 
the outer and inner boundaries. 

4 The composition of the ash has negligible effect on the thermal structure 
of the fuel layer. 
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Fig. 2. — Burst energy and ignition column density as a function of the 
outward flux from the neutron star crust. The boxes mark the regions where 
the theoretically predicted burst energy matches the observed energy for 4U 
1820-30 (green color) and intermediate long bursts (magenta color). The line 
styles are the same as in Fig. [T] 

precise value of Qi, is sensitive to the poorly constrained neu- 
trino emissivity in the core. In this work, we take Qb as a free 
parameter, with the constraints pr ovided by parameter st udies 
of the neutrino cooling efficiency (ICumming et al.ll2006h . We 
iterate the integration process by adjusting the surface flux to 
match the inner bo undary condition . This method is similar 
to that described in lCumm ing (2003|). 

We neglect compressional heating ~ CpT a 0.02 T% 
MeV/nucleon, which is negligible compared to crustal 
heating. Also, for simplicity, we ignore the moder- 
ate energy release from He burning during accumula- 
tion (~ 0.03 (Ay/0.05) MeV/nucleon with Ay being the 
mass fraction of burned He prior to ignitionfl There- 
fore, the thermal structure of the accumulated enve- 
lope is uniquely determined by the outward flux Fb = 
lO^ergcnrV 1 (Q b /0A MeV)(m/10 4 g cm" 2 s" 1 ), a prod- 
uct of Qb and m. For each Fb, we can determine the 
ignition column density y lgn by solving the thermal equa- 
tions. Once yig n is determined, we estimate the energy 
of the burst by E b = 47r^ 2 y ign g nuc /(l + z) = 1.53 x 
10 41 ergs |jign/(10 10 g cm~ 2 )]/(l + z). Here we assume that 
all of the accumulated fuel is burned in the burst and that the 
energy release from complete He burning up to the iron group 
is Gnuc = 1 .6 MeV/nucleon. 

In Figure [2] we show the ignition column density and the 
energy of the burst as a function of Fb- The new 3a rate of 
OKK09 significantly reduces both the ignition column den- 
sity and the burst energy, especially for accreting neutron stars 
with low crustal heat fluxes (at fixed Qb, this corresponds 
to low mass accretion rates). The burst energy is reduced 
by a factor of 10 (100) at F h = 10 21 (10 20 )ergcm" 2 s ' 
~ 7 x 10 39 ergs (7 x 10 40 ergs). 



to 



5 If more than 10% of the accumulated He is burned prior to ignition, 
the He burning energy release could play a role in heating the accumulated 
matter. This would lead to a smaller He ignition column density. 



4. APPLICATION TO ULTRA-COMPACT X-RAY BINARIES 

Ultra-compact X-ray binaries (UCXBs), typically charac- 
terized by a low optical to X-ray flux ratio, have very short 
orbital periods of less than ~ 1 hr. The small orbital separa- 
tion implies that the companion is an evolved star that has lost 
its envelope and is quite hydrogen poor. 

4U 1820-30, located in the glob ular cluster NGC 6624, has 
an orbital period of 1 1 .4 minutes (IStella et alJI 19871) and ma y 
have a helium white dwarf secondary (iRappaport et al.lfl987h . 
It shows frequent type I X-ray bursts with total burst energy of 
« 2-4x 10 ergs, occurring at a local mass accretion rate m « 

0. 2mEdd = 3 x 10 4 g cm" 2 s" 1 . The ratio of persi stent fluence 
to burst fluence is a a 120 (lHaberlet all 1 19871) . suggesting 
burning of He-rich material. The crustal heating Qb at this 
mass accretion rate varies in a range of 0. 1 - 0.4 MeV/nucleon 
for a wide range of core neutrino emissiviti es based on the 
self-co nsistent theoretical ignition models of ICumming et alj 
(2006). Therefore, F b * 3 - 10 X 10 21 erg cm" 2 s" 1 . The 
flux range together with the range of observed burst energies 
constrains the region where the theoretical predictions match 
observations. This region is marked by a green-colored box in 
Fig. 12 All 3a reaction rates can explain the observed energies, 
but the Fy05 and FL87 rates require high crustal heat fluxes. 

There are three UCXBs that have shown intermediate long 
bursts with durations of 10-40 min and burst energies 
of « 0.5 - 2 x 10 41 ergs, occurring at local mass accre- 
tion rates of m « 0.01 WEdd = 1.5 X 10 3 g cm" 2 s _1 . 
Such intermediate long bursts have durations and energies 
intermediate between regular type I X-ray bursts and super- 
bursts with durations longer than an hour. They may be due 
to explosive burning of a large pile of He slowly accumu- 
lated prior to the burst (ICumming et alJ l2006: in 't Zand et alj 
120051). Sources that have shown such bursts are SLX 1737- 
282 ( in'tZandet^l2^p lM angaetal.ll2008l). SLX 1735- 
269_ Molkov et all 120051) and 2S 0918-549 lrn 't Zand et alj 
120051) . SLX 1737-282 has to date shown exclusively interme- 
diate long bursts with a recurrence time of ~ 90 days. It is 
also the only source that has shown more than one intermedi- 
ate long burst. 

The uncertainty in the outward flux from crustal heating is 
dominated by the uncertainty of the mass accretion rate for 
intermediate long bursts. Given an uncertainty of a factor 
of 2 in m, i.e., m ta 750 - 3000 g cm" 2 s" 1 , corresponding 
to a variation of Qi, in the range of 0.9 - 1.4 Mev/nucleon 
and 0.4 — 1.2 MeV/nucl eon at the minimum a nd maximum 
mass accretion rate (see ICumming et all 120061) . the relevant 
range for the outward flux Fb is 0.7 - 4 x 10 21 erg cm" 2 s" 1 . 
The magenta-colored box in Fig. [2] marks the constraint on 
the properties of intermediate long bursts. In the constrained 
range of outward flux, the new 3a rate predicts burst energies 
< 10 40 ergs, which is much lower than the observed energies 
of intermediate long bursts. 

At very low mass accretion rates (m < 0.005 /«EddX al- 
most all of the deep crustal heating reaches the surface, 

1. e., Qb - 1.4 MeV/nucleon. This leads to F b = 2 x 
10 20 erg cm" 2 s" 1 (m/0.001 mEdd)- In this regime, the He ig- 
nition column density and predicted burst energy differ by or- 
ders of magnitude between the OKK09 and Fy05/FL87 3a 
rates. For instance, at m = 0.001 mEdd, the Fy05 rate pre- 
dicts a burst of 5 x 10 42 ergs (comparable to the energy of 
a superburst) while the OKK09 rate predicts a burst energy 
of only ~ 4 x 10 40 ergs. He ignition at a deep column den- 
sity of yign ~ 10" gem -2 could be a possible cause of su- 
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perbursts (see also iKuromizu et all 120031) and has been pro- 
posed as a possible explanatio n for a superburst from 4U 
06 1 4+09 1 ( iKuulkers et al J 120 1 Oh . Therefore, any superburst 
detected at such low accretion rates can be used to constrain 
the 3a reaction rate at low temperature and high densities. 
However, this kind of superburst is likely to be very rare since 
sources with such low accretion rates might not b e able to 
sustain continuous accretion (e.g. Lasota et al. 2008). Even if 
they are persistently accreting, it takes nearly ~ 100 years to 
accumulate enough fuel to trigger a superburst. 

5. DISCUSSIONS AND CONCLUSIONS 

We have investigated the He ignition conditions on accret- 
ing neutron stars b y considering a new 3a reaction rate of 
lOgata et ail (120091) that is higher by orders of magnitudes 
at low temperatures than suggested by previous work. We 
demonstrate that the He ignition condition on neutron stars 
is sensitive to the choice of the 3a reaction rate. The new 
rate leads to considerably smaller ignition column density 
and, hence, lower burst energy. For m > 0.001 niEdd, the 
OKK09 rate predicts a maximum He ignition column den- 
sity of ~ 3 x 10 9 g cirT 2 , corresponding to a burst energy of 
~ 4 x 10 40 ergs. We have compared theoretical predictions 
of burst properties obtained with multiple 3a rates to obser- 
vations of intermediate long bursts in UCXBs, using the as- 
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sumption that these s ystems are pure He accretors. We find 
that the new 3a rate of lQgata et alj J2009I) is inconsistent with 
observations and cannot explain the observed energy of such 
intermediate long bursts while calculations with the previ- 
ous rates (e.g. iFynbo et al.ll2005HFushiki & Lambll987l) yield 
predictions in general agreement with burst observations. 

Some binary evolutionary models predict that UCXBs 
might have a small amount of hydrogen (X ~ 0.1) in the 
accreted matter (e.g. Podsiadlowski et al. 2002). Presence of 
hydrogen provides additional heating from hydrogen burning 
during accumulation which would lead to even smaller igni- 
ti on column d e nsities (ICummin gl l2003l) . Hence, the 3 a rate 
of lQgata et al.l J2009I) would have an even harder time to ex- 
plain intermediate long bursts in systems that are not pure He 
accretors. 



It is a pleasure to thank Edward Brown for helpful com- 
ments which improved the paper. We thank Evan O'Connor, 
Hendrik Schatz and James Truran for useful discussions and 
comments. We thank the anonymous referee for helpful sug- 
gestions. We also thank Kazuyuki Ogata for providing us with 
his triple-a reaction rates. This work is partially supported 
through NSF grants AST-0855535 and OCI-0905046. 



Kuulkers, E., et al. 2010, A&A, 514, A65+ 
Lasota, J., Dubus, G, & Kruk, K. 2008, A&A, 486, 523 
Molkov, S., Revnivtsev, M., Lutovinov, A., & Sunyaev, R. 2005, A&A, 434, 
1069 

Narayan, R., & Heyl, J. S. 2003, ApJ, 599, 419 

Ogata, K, Kan, M., & Kamimura, M. 2009, Progress of Theoretical Physics, 
122, 1055 

Peng, F, Brown, E. F, & Truran, J. W. 2007, ApJ, 654, 1022 
Podsiadlowski, P., Rappaport, S., & Pfahl, E. D. 2002, ApJ, 565, 1107 
Potekhin, A. Y., & Chabrier, G. 2000, Phys. Rev. E, 62, 8554 
Potekhin, A. Y., Chabrier, G, Chugunov, A. I., Dewitt, H. E., & Rogers, 

F. J. 2009a, Phys. Rev. E, 80, 047401 
Potekhin, A. Y., Chabrier, G, & Rogers, F. J. 2009b, Phys. Rev. E, 79, 

016411 

Rappaport, S., Ma, C. P., Joss, P. C, & Nelson, L. A. 1987, ApJ, 322, 842 
Saruwatari, M., & Hashimoto, M. 2010, AIP Conference Proceedings, 1269, 
412 

Sato, K. 1979, Prog. Theor. Physics, 62, 957 

Stella, L., White, N. E., & Priedhorsky, W. 1987, ApJ, 312, L17 



